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ABSTRACT 

The Extremely Large Telescopes currently under construction have a collecting area that is an order 
of magnitude larger than the present largest optical telescopes. For seeing-limited observations the 
performance will scale as the collecting area but, with the successful use of adaptive optics, for 
many applications it will scale as (where D is the diameter of the primary mirror). Central to 
the success of the ELTs, therefore, is the successful use of multi-conjugate adaptive optics (MCAO) 
that applies a high degree correction over a field of view larger than the few arcseconds that limits 
classical adaptive optics systems. In this letter, we report on the analysis of crowded field images 
taken on the central region of the Galactic globular cluster NGC 1851 in Kg band using GeMS at the 
Gemini South telescope, the only science-grade MCAO system in operation. We use this cluster as a 
benchmark to verify the ability to achieve precise near-infrared photometry by presenting the deepest 
Kg photometry in crowded fields ever obtained from the ground. We construct a colour-magnitude 
diagram in combination with the F606W band from HST/ACS. As well as detecting the “knee” in 
the lower main sequence at Kg :^20.5, we also detect the double subgiant branch of NGC 1851, that 
demonstrates the high photometric accuracy of GeMS in crowded fields. 

Keywords: globular clusters: individual (NGC 1851) — instrumentation: adaptive optics — tech¬ 
niques: photometric 


1. INTRODUCTION 

In classical adaptive optics (AO), the aberrated wave- 
front of a guide star is measured by a wavefront sensor 
and corrected by a deformable mirror (DM) to deliver a 
diffraction-limited image inside a field of view only sev¬ 
eral arcseconds wide, outside of which the point spread 
function (PSF) rapidly widens and deforms. This scale 
is set by the limiting angle inside of which the turbu¬ 
lence is correlated. A larger and more uniform field of 
view can be obtained with multi-conjugate adaptive op- 
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tics (MCAO), by using an asterism of guide stars with 
multiple DMs optically conjugated at different heights of 
the atmosphere. Typically, AO works best in the near- 
infrared (NIR) because the spatial and temporal coher¬ 
ence of the aberrated wavefront is larger at longer wave¬ 
lengths, resulti ng in better corrections (f or a comprehen¬ 


sive review, see 


Davies & Kasperl (|2012|)). 


The Gemini Multi-Conjugate Adaptive 0 

ptics System 

(GeMS) (Rigaut et al.| 

2014 

Neichel et al.| 

|2014) on the 

Gemini South telescope on C 

Jerro Fachon in Chile is the 


first facility-class MCAO instrument in operation, al¬ 
though the concept has previously been proven with the 
Mnlti-C onjngate Adaptive Op tics Demonstrator (MAD) 
on VLT (Marchetti et al. 2007|. GeMS uses two DMs and 
five sodium laser guide stars (LGS) to deliver a corrected 
field of view of 83" x 83" to the Gemini South Adaptive 
Optics Imager (GSAOI) ( McGregor et al.|2QQ4 i Carrasco 
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Figure 1. Left: Mosaic image of NGC 1851 taken with GeMS/GSAOI with a field of view of 83^^ x 83^^ The positions of the five sodium 

are shown witJi orange markers and indicated by red circles are the three lN(Jbs (clockwise trom the iett, rriR = l'6.'6, lAA, 1'6. (). 

tCiQtit: a zoom-in ot tJie region in the lett panel contained witJiin the blue rectangle observed witii UbAUl in Ks band (toy) and the same 

mTEe^pncaO^omparaBIe^^^^^^™ — 


et al.||2Q12 |, a mosaic of 2 x 2 HAWAII-2RG chips with 


a pixel scale of 0.01962" px'^ 

Galactic globular clusters (GGGs) are an ideal target for 
MCAO systems. From an instrumentation perspective, 
their moderately large expanse on the sky and their high 
degree of crowding means that they are well-suited to 
MCAO systems that can potentially resolve faint stars in 
their central regions, while still possessing a large num¬ 
ber of bright red natural guide stars (NGS) that are 
needed for the determination of the tip-tilt and defo- 
cns a berrations (jRigaut & Gendron||1992| |Herriot et al. 

far photome- 

try of GG Gs has been a prerogative of the Hubble Space 
Telescope (Milone et al. 2012 2014), in the future AO- 
assisted EL'l's could get an edge over space telescopes in 
the NIR due to their much smaller diffraction limit. 
From an astronomical perspective, the last decade has 
witnessed a revolution in our understanding of GCs. The 
current spectroscopic evidence indicates that all GGGs 
show star-to-star variations of C and N abundances to¬ 
gether with well defined 0-Na and Mg-Al anticorrela¬ 
tions in both evolved (RGB, HB) and unevolved stars. 
The fraction of GGGs showing also multiple sequences is 
limited and probably correlated with the total mass of 
the cluster (jGratton et al. 2012). This confirmation of 
multiple stellar populations is m contrast to the classi¬ 
cal view of a coeval cluster with homogeneous chemical 
abundances. In the case of NGG 1851, a doubl e subgiant 


bra nch (S GB) has famously been observed by Milone et 


al. (2008) using HST/AGS images. This feature has been 


interpreted as evidence of a second generation of s tar for¬ 
mation ^ 1 G yr older than the first; alternatively, [Gassisi 


et al. (|2008) suggested that the double SGB can be ex- 


plained as two coeval populations, where one is a factor 
two more enhanced in CNO than the other. 

Here, we present first results from a Gemini/GSAOI cam¬ 
paign targeting Galactic GGs using GeMS. In particular, 
we present observations of NGC 1851 taken in the Kg 


band that represent the deepest Kg photometry of a GC 
taken from the ground. We use this GC as a bench¬ 
mark to demonstrate the precision of the MCAO data 
by producing a CMD of the stellar population with an 
HST/ACS optical band. 

2. DATA 

During the GeMS/GSAOI Science Verification phase 
in Fall 2012, we obtained images of NGC 1851 in the J 
(A = 1.235 /im) and Kg (A = 2.159 /im) bands. Subse¬ 
quent observing runs targeted five other globular clusters 
in the same bands. For this first paper in this series, we 
describe analysis of the Kg band data for NGC 1851; a 
full analysis of the J and Kg band discussing the prop¬ 
erties of NGC 1851 will be published in a forthcoming 
article, and the other GCs will be discussed in subse¬ 
quent papers. The observation with MCAO of globular 
clusters in NIR will prove particularly useful to study 
low mass stars, giving information on the total cluster 
mass, the lower end of the mass function and the colour- 
temperature relation for cold atmospheres. We will also 
discuss the use of these data for astrometric analyses. 
The median seeing of our observations was recorded by 
RoboDIMM on Cerro Tololo to be of 0.77" and the in¬ 
strument delivered in Kg band an average corrected full 
width at half maximum (FWHM) of 0.09", comparable 
to the 0.06" FWHM of the diffraction limit case. The 
integration consisted of 12 dithered sub-exposures of 160 
s each, and the time from the first to last science image 
was 75 minutes. Two short exposures of 21 s and 90 
s were observed to obtain photometry for the brighter 
sources that were saturated in our longer exposures. 
Image reduction has been performed using the Gemini 
GSAOIIRAF package that contains a bad pixel map and 
appropriate non-linearity correction. A super-flat was 
made that combined the twilight and dome flat fields, 
utilising the uniform spatial illumination of the former 
and the high signal-to-noise ratio of the latter. The low 
dark current of the H2RG detectors means we did not 
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Figure 2. Left: (Kg, F606W-Ks) colour-magnitude diagram of NGC 1851; the detail of the double SGB is shown in the inset. Right: same 
as the left panel with average photometric (random) uncertainties indicated. Overlaid is the fiducial line with the approximate locations 
of the main sequence turnoff and main sequence knee highlighted by red dots. 


need to apply a dark frame. Local sky estimates were 
used for the point source photometry (see below), and so 
we did not need to apply a sky frame either. However, it 
is worth noting that there is a substantial number of bad 
pixels in each chip, and that one of the chips has a ^ 0.2 
mag higher efficiency than the other. Figure shows an 

example of a processed image. __ 

The DAOPHOT suite of progr ams ( Stetson|1987| [Stetson 


& Harris 1988 Stetson|1994 ) has been used to perform 


the photometric analysis on the individual frames , hav¬ 


ing already been proven to work with MCAO data ( Bono 


et al. |2Q1Q[ IFiorentino et "ah 2011 |Davidge||2014 


Sara- 


cino et al.||2015| ). To remove cosmic rays and bad pixels 


from our catalogue, we keep only objects that are found 
within a radius of one pixel in at least three dithered 
exposures out of twelve. Particular care has been paid 
to the PSF fitting routine due to the expected shape 
of the MCAO PSF. AO concentrates the light into the 
core of the PSF, but only up to a certain scale given by 
the density of actuators i n the DM, beyond wh ich the 
PSF remains uncorrected ( [Jolissaint et al.pQQG ). Thus, 
even though the fraction of light in the core of the PSF 
is substantially greater than in the seeing-limited case, 
the size of the halo relative to the core is much larger 
in the AO case, and the adopted PSF model should re¬ 
flect this. In addition, even although MCAO provides a 
more uniform correction over a larger field than classi¬ 
cal AO, some significant spatial and temporal variation 
in the PSF structure may still be expected. Therefore, 
each exposure has ^ 300 manually selected bright and 
isolated stars used to model the PSF as an analytic pro¬ 
file plus a spatially-varying look-up table, where each 
pixel of the PSF is allowed to spatially vary as a cubic 
function, independent of every other pixel of the PSF. 
While we conduct our analysis independently on every 


exposure of every chip, we also perform an analysis on 
the stacked median image of each chip in order to detect 
and photometer the faintest stars. We ensure a high fi¬ 
delity for our stellar catalogue by cross-matching every 
detection to the catalogue of NGC 18 51 from the ACS 
Surve y of Galactic Globular Clusters (Sarajedini et al. 
2007|). Noise features and artefacts that are incorrectly 


identified as stellar objects are removed by this process, 
while real stars are kept. The full HST/ACS catalogue 
is considerably deeper than the Kg catalogue (by ^4.5 
magnitudes), and so this does not limit the depth of the 
Kg band data presented. 

Each frame of each chip is calibrated independently to 
the Kg band magnitude scale of the 2MASS Point Source 
Catalog through an intermediate catalogue generated by 
one of the authors (PBS) from 256 seeing-limited archival 
images obtained with CTIO/NEWFIRM. These latter 
exposures are both wide and shallow enough to contain a 
large number (^ 4000) of stars present also in the 2MASS 
catalogue. These were used to calibrate the magnitudes 
of fainter stars in the same images, which could then 
be used to determine the photometric zero points of the 
GeMS images. 

We note that the extended halos of the MCAO PSE, 
combined with the spatial variation in the shape of the 
PSE, can lead to a photometric zero-point that changes 
between sub-exposures. This will be discussed and ex¬ 
plored in more detail in Turri et al. (2015, in prepara¬ 
tion). The data presented in this letter have been cor¬ 
rected for this variable zero-point, but no color correction 
is applied, however we estimate it to be small. Note that 
this makes no difference to the analysis presented herein. 
Eor each star, the final magnitude was calculated by com¬ 
bining the magnitud es measured in each exposure using 
artificial skepticism (Stetson 1989 chap. 3). The red- 
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Figure 3. Left: As Figure but enlarged in the region of the double SGB. Overlaid is the selection box orientated parallel with the linear 
segment of the SGB. RighLnistogrSiin of the distribution of offsets relative to the median of the left panel box. Overlaid is the best-fit 
double Gaussian to the bimodal distribution. 


dening correction i n the d i rection of NGC 185 1 is small 
{E{B -V)= 0.02, 
not been applied. 


Harris 


( 1996, 2010 editiolil )) and has 


3. RESULTS 

Figure 1^ shows the (Kg, F606W-Ks) CMD for NGC 
1851 obtained by our analysis, and contains ^ 16000 
stars. The error bars are a combination of random uncer¬ 
tainty (photon and readout noise), PSF fitting residuals 
and scatter of the measurem ents between exposures . The 
F606W data are taken from Sarajedini et al. (2007), and 
the Kg data have additionally had quality control criteria 
in sharpness and applied. In addition, a coarse proper 
motion cut has been applied to reject foreground stars by 
restricting the apparent difference in the positions of the 
stars between the two epochs to 8 mas (the HST/ACS 
data were observed in 2006). More details and astromet¬ 
ric analysis will be presented in Turri et al. (2015, in 
preparation). 

Several notable features are visible in Figure Firstly, 
the extended horizontal branch is clearly visible and the 
red HB is particularly well defined. The asymptotic gi¬ 
ant branch can be seen departing from the red HB just 
prior t o the point where our data are saturated. In ad¬ 


dition, Saviane et al.| ( 1998 ) have observed the red gi¬ 
ant branch bump in NGC 1851 at a magnitude similar 
to the red HB. This feature occurs when the H-burning 
shell moving outwards approaches the chemical discon¬ 
tinuity caused by the first dredge-up of the convective 
envelope, and it can be observed in Figure on the RGB 
at Ks=13.8. 

The depth of the Kg is of particular note. Figure 
reaches ^3.5 magnitudes deeper than the main sequence 
turnoff with a magnitude error of 0.15 at Kg=22 using 
a total exposure of 1920 s. This depth is impossible to 
achieve for a crowded target like NGC 1851 without the 
us e of MG AO. A s imilar depth in Kg has been obtained 


by Monelli et al. (2015) in M 15 using FLAG on LBT 
but 3' from the center of the cluster. 

A bend in the main sequence is visible towards the limits 
of our data, and starts at Kg 20.5. This feature, known 
as the main sequence knee (MSK) is a typical feature in 
NIR colour-magnitude diagrams of GCs and is caused by 


cold stellar atmospheres ( 

Bono et al.|2010| 

Monelli et al. 

2015). It has been proposed that the rel 
between the MSK and the main sequence 

ative distance 
turnoff can be 


used as a dista nce and reddening inde pendent cluster age 
determination (Calamida et al.|2009), being the position 


of the former point independent of age. We will revisit 
this claim in future papers in this series. 

Of particular interest, the subgiant branch in Figure 
is clearly split into two seque nces. This fea t ure w as first 
identified for NGC 1851 by |Milone et al.| (|2008|) using 
HST/ACS data. Previous observations of the core of a 
GC with double SGBs from the ground have revealed 
at most a broadened SGB with a bimodal distribution 


(specifically, the ex cellent study of NGC 6656 by Li- 
bralato et al. ( 2Q14| )); only observations of the uncrowded 


outskirts have been a ble to discern multiple SGBs from . 
seeing-limited images ([Milone et al.|2Q Q9 Han et al. 2009 

d. T'he fact tha 


Cummings et al.||2014j). T'he fact that this feature is vis- 

ible in the centre of JNGC 1851 using the Kg band data 
from GeMS is a verification of our small photometric 
uncertainties which demonstrates significant promise for 
the successful exploitation of MCAO for precision pho¬ 
tometry of crowded fields in the ELT era. 

The right panel of Figure shows a histogram of the 
offsets of the stars contained within the box in the left 
panel. Two separate peaks are clearly identifiable and 
these correspond to th e faint and b r ight S GB (fSGB, 
bSGB) as described by |Milone et al.| (|2008|). We model 
this distribution as two Gaussian curves that are over¬ 
laid on the histogram. By integrating these curves, we 
find Nfso B/NhsoB = 0 - 47, en tirely consistent with the 
findings of Milone et al.l (|2009|). 


4. GONGLUSION 


We present first results from a study of the core of NGC 
1851 using the GeMS MCAO system on Gemini South. 
We demonstrate the ability of this system to achieve pre¬ 
cise near-infrared photometry of crowded regions by pre¬ 
senting the deepest Kg observations ever obtained from 
the ground. In the colour-magnitude diagram created 
with the addition of the HST F606W band, we detect 
the “knee” in the lower main sequence and the double 
subgiant branch of the globular cluster. A subsequent 
paper will present a more detailed analysis of NGC 1851 
using these data, as well as an analysis of the photometric 
and astrometric performance and stability of the GeMS 
system. Here, we have used this cluster as a benchmark 
of the quality of the MCAO photometry that has been 
obtained. This work is part of a larger project with 
GeMS to obtain J, Kg photometry for several GCs in 
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the southern hemisphere some of which have NIK da ta 
from HST/WFC3, e.g. NGC 2808 (|Milone et al.||2012D. 
MCAO systems are an integral part of the two largest 
ELT projects as they tackle science programs that re¬ 
quire an extended field of view, with NFIRA OS as a first 
light system for the Thirty Meter Telescop e dBoyer et al 
2014 ), and MAORY on the European-ELT (Diolaiti et al 


2014). As the only science-grade MCAO system, GelVlS 


IS an important new facility instrument that can serve as 
a scientific testbed for developing a practical understand¬ 
ing of the challenges of obtaining accurate photometry 
and astrometry using such a system. This pilot study 
implies that high precision photometry can be achieved 
from the ground using MCAO, and can inform the de¬ 
velopment of appropriate data processing techniques and 
observing strategies to ensure the ELTs deliver their full 
scientific promise over extended fields of view. 


Based on observations obtained at the Gemini Obser¬ 
vatory, which is operated by the Association of Univer¬ 
sities for Research in Astronomy, Inc., under a cooper¬ 
ative agreement with the NSE on behalf of the Gemini 
partnership: the National Science Foundation (United 
States), the National Research Gouncil (Ganada), CON- 
IGYT (Chile), the Australian Research Council (Aus¬ 
tralia), Ministerio da Ciencia, Tecnologia e Inovacao 
(Brazil) and Ministerio de Ciencia, Tecnologfa e Inno- 
vacion Productiva (Argentina). Acquired through the 
Gemini Science Archive and processed using the Gemini 
IRAF package. 

G. Fiorentino has been supported by the FIRB 2013 
(MIUR grant RBFR13J716). 

We appreciate the prompt support received from the 
GeMS staff for the image reduction. We are grateful to 
the adaptive optics group at NRC Herzberg for feedback 
on the system performance and to INAF - Osservatorio 
Astronomico di Bologna for insights on stellar popula¬ 
tions in globular clusters. 

Facilities: GeminiiSouth (GeMS, GSAOI). 
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